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The problem of dissipative motion of solitons in the case of tetramethyl ammonium manganese 
chloride (TMMC) is studied as a function of the external magnetic field and the temperature. Two 
specific situations are analyzed separately; the first, above the transition temperature Tjv, in which 
the classical motion of the spin degree of freedom is described by a sine-Gordon (SG) equation of 
motion and, the second, below Tjv, in which the system is described by a double sine-Gordon (2SG) 
equation of motion. The existence of a dissipative regime for the soliton motion and its influence 
on the dynamical structure factor -which might be experimentally detected- are reported. 
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I. INTRODUCTION 



In the past few decades it has become well-established 
that the physical properties of some magnetic materi- 
als, TMMC, CsNiF 3 (caesium nickel fluride) and CuCl 2 
2NC5H5 (dicloro-bis-piridine copper II), for instance, 
have essentially one dimensional character above their 
transition temperature. u In those kind of materials the 
distance between magnetic ions along a given direction 
(magnetic chain direction) is shorter than in the other di- 
rections. In such an arrangement the intrachain coupling 
constant is typically more than two orders of magnitude 
stronger than the interchain coupling constant. There- 
fore, the system can be considered as a set of weakly 
interacting magnetic chains . Due to the relative sim- 
plicity of obtaining solitonic or solitary-wave solutions in 
ID systemsB, these quasi-one-dimensional magnets turn 
out to be the paradigm for the study of the influence of 
the non-liner modes (solitons) in the dynamical proper- 
ties of such systems at finite temperatures.™ Although 
all real magnetic materials investigated are not perfectly 
one dimensional, the assumption of the ID behavior is 
shown to be in good agreement with the experimental 
results (see Ref.l and the references therein). 

In magnetic materials solitons or solitary- waves can be 
regarded as 'kinks' or 'twists' in the spin space moving 
with constant speed and carrying a constant topological 
charge defined by the values of the spin variables at infin- 
ity. For low enough temperatures, when the linear modes 
(spin-waves) are not excited, the magnetic system can 
be represented in fixst approximation by a gas.-pf non- 
interacting solitons.B Using this idea, MikeskaoQ calcu- 
lated the soliton contribution to the dynamical structure 
factor of the classical one-dimensional magnets. From 
both works we learn that the assumption of ballistic mo- 
tion for solitons is the origin of the 'central peak' behavior 
observed in neutron scattering experiments. A different 
situation could be found from the quantum field theory 



point of view when the temperature is raised. In this case 
the spin-wave (SW) modes are excited, therefore not all 
of the degrees of freedom of the system contribute to the 
soliton formation and a residual interaction (which cou- 
ples the center of mass of the soliton to the spin-wave 
modes) shows up. In practice, the specific form of this 
kind of interaction is obtained via the collective coordi- 
nate methodtl in the quantization process! of the classical 
hamiltonian. 

The soliton-SW coupling may result in a dissipative 
regime to the soliton motion depending on the form of 
the potential generated by the presence of the non-linear 
excitation.EEl As it is known, the equation of motion for 
the spin variable in the TMMC below and above Tjv are 
a 2SG and a SG equation^ respectively. This fact makes 
the TMMC a suitable probe to investigate the appear- 
ance of a dissipative regime in the soliton motion pro- 
vided that, below and above a certain Neel transition 
temperature T/v, the classical equation of motion for the 
spin variables are substantially different. The main pur- 
pose of this work will be therefore the analysis of the mag- 
netic soliton motion above and below T/v and the possible 
influence of the dissipative regime, found for T < T/v, on 
the solitonic contribution to the dynamical structure fac- 
tor. In doing that, we will use the method developed in 
Ref. 9 for the analysis of the dissipative dynamics of 
solitons. As pointed out before, this -formalism is based 
on the collective coordinate methodtl and allows us to 
transform the original hamiltonian of the spin degree of 
freedom into one of a particle (the soliton) coupled to an 
infinite set of linear- modes (SW). 

Starting from the interacting soliton-SW hamiltonian 
it is possible to obtain a Brownian-like equation of motion 
for the soliton center of mass via the Feynman- Vernon 
formalism. This effective equation of motion is writ- 
ten in terms of a damping constant that depends on the 
phase shifts of the scattering problem that emerges by 
the presence of the non-linear excitation coupled to the 
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SW. Therefore, the analysis of the scattering properties 
of the 2SG potential (for T < T N ) and the SG potential 
(for T > Tjv) allows us to calculate the mobility as a 
function of the temperature and the external magnetic 
field. As it will be shown, above Tjv the SG solitons have 
infinite mobility in agreement with the ballistic motion 
used to understancLthe neutron scattering experiment for 
H/T > lOfcOe/idil On the other hand, for low temper- 
atures or low magnetic field, when the spin equation of 
motion for TMMC have a 2SG form, the soliton mobility 
is finite, changing the form of the dynamical structure 
factor considerably. 

To begin with, in Sec. II we review the models cur- 
rently applied to the spin dynamics of the TMMC com- 
pound, above and below its transition temperature, and 
also the corresponding classical equations of motion. In 
Sec. Ill we summarize the obtainment of the quantized 
soliton-SW hamiltonian and the damping parameter of 
the soliton Brownian motion is calculated as a function 
of the temperature and the external magnetic field. Sec. 
IV is devoted to the study of the influence of the soliton 
damped motion on the dynamical structure factor and, 
finally, our conclusions are presented in Sec. V. 



II. THE MODEL FOR TMMC 

The antiferromagnet TMMC has extensively been 
studied from the theoretical and experimental points of 
view. The hamiltonian describing the interacting 3D ar- 
ray of classical spins in this material can be written as 



hjSi'j, 



(1) 



ijti' j 



where 



H i = E {J\\S^S hk+1 + A(Sl k ) 2 - gnBBSf ik } . (2) 

k 

The hamiltonian Hj describes the nearest neighbour in- 
trachain interaction between spins with an easy plane 
anisotropy (A > 0) placed in an external magnetic field 
(B) in the x direction. The spins will be treated as classi- 
cal vectors of lengh S and the constants J|| and J± , both 
positive, correspond to the antiferromagnetic and ferro- 
magnetic exchange coupling constants, respectively. The 
second term in the r.h.s. of (|l|) represents an interchain 
interaction between the spins, completing the description 
of the 3D spin arrangement. Finally, the following val- 
ues of material parameters will be used: Jii = 13. AK, 
S = 5/2, A/J|| = 0.01 - 0.02, Jj_/J|| = 1.5 • 10~ 5 and 
g = 2.01. 

In order to start the classical description of the spin 
dynamics it is convenient to look at two main different 
situations, namely, temperatures below and above the 
transition temperature. For temperatures below Tjv the 
system described by d|) displays a long range magnetic 



order, therefore the staggered spontaneous magnetization 
is not zero and the system can be described in the mean 
field approximation as a set of non-interacting antiferro- 
magnetic chains with an additional spontaneous magne- 
tization in the y directiontx Explicitly, 

H = E { J n^^+i + A (s?) 2 - g^sBS't 

i 

-g^ B B^(-iysy}, (3) 



where 



B 



MF 



r,J ± ((-iySf)/g^ 



(4) 



and rj accounts for the presence of neighbouring chains 
in the model. In the specific case of TMMC, rj = 6. 
The intrachain mean field B^ F is usually replaced by 
its saturation value Bf_ ss 22.30e which results from the 
substitution of s[ v) in (@) by its maximum value. 

At this point, we can carry on the classical description 
of the spin dynamics. In order to do that it is convenient 
to change the spin variables to the following form 



S e , = ±S [sin(6 ± 6>) cos($ ± (p), 

sin(9 ± 6) sin($ ± <p), cos(6 ± 6)] 



(5) 



where e and o stands for even and odd sites within a 
chain. 

Using the representation (|5|) a ^-dependent part of the 
hamiltonian (g) can be obtained (see Ref. 3 for details). 
Explicitly, 



-JnS 2 



dz 



;(d t <i>) 2 + (9 Z $) 2 



1 



-b sin 2 $ - 2b± sin$ 



where 



c 2 = 4 
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J \\ 



9^bB 
JnS ' 



6j = 



9^bB 



MF 
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(6) 
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It should be stressed that the hamiltonian (^|) is an ap- 
proximated description of the real TMMC system. To re- 
produce the experimental results, magnon-mass and soli- 
tonic energy, for instance, quantum effects and the out 
of plane component of the magnetization must be taken 
into account £3 To go on with the classical description of 
the ^-dependent part of the original hamiltonian ([!]) the 
equation of motion associated to (pi), 



1 b 2 
—d t t® = d zz § sin 2$ - b± sin $ 



(8) 



has to be solved. Equation is not completely intc- 
grable, however, it has solitonic solutions in the form of 
27r-kinks(antikinks) moving with velocity u. Explicitly, 



cos$ = ±2- 



1 + a sinh y 



sinhy, 



(9) 
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sin <j> = 1 — 



where 



6 2 /4' 



1 + a sinh y 



y = {z- z{t)) 




(10) 



(ii) 



and the position of the soliton center of mass z(t) is given 
by 



z{t) = zq + ut. 



(12) 



On the other hand, for temperatures above T/y the 
value of the b± is very small. In fact, in this situation b± 
can be set equal to zero and, a well-known solitonic solu- 
tion for equation (^) can be found: the 7r-kink(antikink) 
solution for the SG equation 



sin$ s (z,i) = ±tanh [(1 - u 2 /c 2 y 1/2 (z - z{t))b/2 



(13) 



As it can be seen, the model for TMMC in the con- 
tinuum approximation leads us to different kinds of soli- 
tonic solutions depending on the temperature. A 2SG 
soliton solution given by @ and ( [Tof ) for T < and, 
a SG solution (|l3|) for temperatures above TV. As it 
was already mentioned, from the classical point of view, 
these soliton solutions will move with constant velocity 
througout the sample. However, looking at the soliton 
dynamics from the quantum field theory perspective, the 
interaction with the spin waves can transform this bal- 
listic regime into a dissipative one. In the next section 
we shall be aiming at the investigation of the mobility of 
the two types of solitons, below and above TV- 



where the functional integral has the same initial and fi- 
nal configurations and tr means to evaluate it over all 
such configurations. As the functional integral in (|l^) is 
impossible to be evaluated for a potential energy density 
as in (^) we must choose an approximation to do it. 
Since the magnetic moments at the manganese sites in 
the TMMC are large (5/2), the semi-classical limit will 
be chosen as the appropriate one in our case. Within the 
functional integral formalism of quantum mechanics, the 
semi-classical limit is simply the stationary phase method 
applied to ( |l6| ) around the solitonic solutions (|9|) , ([!(]) or 
( |l3| ) in which we are interested. When this is done we 
are left with an eigenvalue problem that reads 



+ U"{$s) \ i> n {z - zq) = k 2 n i> n (z - z ), (17) 



dz 2 



where $ s is denoting the soliton-like solution around 
which we are expanding <I>(z,i) and ip n (z — zq) are the 
spin wave modes in the presence of the soliton. 

Now one can easily show that d$ s /dz is a solution of 
( |l7| ) with k n = 0. The existence of this mode is related 
to the translation invariance of the system and causes 
the divergence of the functional integral in (|l^) in the 
semi-classical limit (Gaussian approximation). The way 
out of this problem is the so-called collective coordinate 
methodu. This method consists basically in expanding 
the field configurations about $ s (z) as 

oo 

$(z,t) = 3> S (Z - Z (t)) Cnlpn (Z - Z (t)) , (18) 

n=l 

but regarding the c-number zq as a position operator. Us- 
ing expansion ([T^|j, the second quantized version of (|^) 
can be written aia 

H = 2W {P ~ £ h9mnb " bm)2 + E m nbtb n . (19) 



III. SOLITON MOBILITY 

The quantum dynamics of our spin system (|^) can be 
analyzed by studying the quantum mechanics of the field 
theory described by the action 



sm = J\\S 2 I I \ 2^(0t<!>! 2 ~<<M>r i n-j,)^/,,/;. 



where 



(14) 



(15) 



To quantize the system described by (|14|) we need to 
evaluate 



[/(<&) = — sin 2$ -I- b± sin 



G(t) = tr / £>$ exp -S[^>] 



(16) 



where f2„ = ck n . 

In the hamiltonian ([l9|), P stands for the momentum 
canonically conjugated to zq, 



2J\\S 2 a f + °° 
M„ = 1L / dzU(*.(z)) 



(20) 



is the soliton masso and the coupling constant g mn is 
given by 



1 



dzipm(z) 



dlpn(z) 

dz 



(21) 



The operators 6 + and b are respectively the creation 
and annihilation operators of the excitations of the mag- 
netic system (magnons) in the presence of the soliton. In 
fact, the term 



E 



(22) 
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can be interpreted as the total linear momentum of the 
magnons of the system and therefore, we are left with a 
problem in which the momentum associated to the soli- 
ton is now coupled to the magnons' momenta. This effec- 
tive model suggests that, as the population of magnons 
is a temperature-dependent quantity, the mobility of the 
soliton will be strongly related to the temperature of the 
system and its dynamics (determined by (|l9|)) will be 
non trivial. 

At this point we are ready to study the mobility of 
the wall because we have been able to map that problem 
into the hamiltonian (|l9|), which on its turn -has been 
recently used to study the mobility of polaronaij, heavy 
particles in lD-environmentali! and skyrmions in 2D elec- 
tronic systemsEj. The main result obtained in those cal- 
culations can be summarized as follows. The damping 
function j(t) (basically the inverse of the mobility) is 
given by 

[n{uj) - n(uj')] cos{u - u')t} , (23) 



where 



1 



(24) 



is the Bose function and, 



S{w,uf) = \9mn\ 2 6(u - n n )S(to' - O m ) (25) 

ma 

is the so-called scattering function. 

In the long time limit j(t) can, to a good approxima- 
tion, be written as 



7 (t) = f(T)S(t), 



(26) 



where 8(t) is the Dirac delta function and j(T) is given 
by 



i 



2ttM s Jo 



dE K{E) 



0E e 



0E 



(eP E - l) 2 ' 



(27) 



In (|27|), TZ(E) is the reflection coefficient of the "poten- 
tial" U"($ s ) in the Schrodinger-like equation (fL7|) . For 
simplicity we will express the reflection coefficient TILE) 
in terms of the even and odd scattering phase shiftstj as 



K{k) = sin 2 {8 e {k)-5°{k)). 



(28) 



At this point we can perform the calculation of the soli- 
ton mobility in TMMC for temperatures above and below 
T N . 



A. Soliton mobility for T <T N 

To calculate the soliton mobility below Tjv we need 
the explicit form of the potential t/"($ s ) involved in (27). 



For the-case of the 2SG soliton this potential can be writ- 
ten asB 



where 



n/^2SG\ 



1 

A 2 



1 - 2sech 2 (- +p)- 2sech 2 (^- - p) 

A A 



Z Z 

+2sech(Y + p)scch(— — p) 
X X 



cosh p 



6 2 /4' 



(29) 



(30) 



The second and third terms in the r.h.s. of (|29j) are 
the potentials of the noninteracting 7r-solitons located at 
z/X = ±p whereas the last term describes the interaction 
of the two 7r-solitons at z/X = ±p respectively. 

For all finite values of A and p, the system is trans- 
lationally invariant and, consequently, the potential (|29| ) 
has a zero energy state that is given by 



ipo oc sech(— + p) + sech(— — p) 



X 



X 



(31) 



which is nothing but the Goldstone mode of the 27T- 
soliton for finite transverse magnetization and finite ex- 
ternal field. 

In order to evaluate the expression (|27j ) for the damp- 
ing constant we need the even and odd phase shifts asso- 
ciated to the potential (^9|). Unfortunately, their analyt- 
ical evaluation is very complicated for all finite values of 
A and p, and in what follows we will only study the sit- 
uation of weak external fields (b± 3> b 2 /2). In this case 
(p<l) the Shrodinger-like equation $LT\) can be written 



(32) 



where 



—^ + V{ Z )\lP n {z) = K 2 n ^ n {z), 



A 2 



A 2 ' 



(33) 



and the potential (|29j) is now reduced to the sum of two 
contributions, one coming from the spontaneous stag- 
gered magnetization and, the other from the presence 
of the weak external field. Explicitly, 



with 



and 



V(z) = V (z) + (^fV 1 (z), 



V {z) = -2sech 2 (|) 
V x {z) = -8tanh 2 (~) scch 2 (£) 



(34) 



(35) 



(36) 



4 



The calculation of the even and odd phase shifts for 
a potential of the form (|3^)-(|3^) is reported in Ref. 10 
and here we will only show the fundamental results of the 
numerical solution of the Schrodinger-like equation (32). 
Fig. 1 and Fig. 2 show the even and odd parity phase 
shifts for different values of the external field. 



Momentum 

FIG. 1. The even phase shift as a function of the momen- 
tum for different values of p. The continuous line correspond 
to p = 0.14, the dotted line to p — 0.31 and the dashed line 
to p = 0.50. 



As it can be seen in Fig. 1 the even phase shift is 7r/2 at 
the origin. This behavior is in complete agreement with 
the existence of an even bound state corresponding to 
the Goldstone mode. On the other hand, the odd phase 
shift 5 o (0) = 7r, indicates the presence of an odd bound 
state. This result was previously obtained by Kivshar et 
aZJ13 in the study of the small-amplitude modes around 
the localized solution of the 2SG equation and shows that 
there is always an odd bound state in this kind of system. 
Therefore, the spectrum of ([54]) is composed by: i) the 
ipo solution ( pil] ) corresponding to the translation mode of 
the soliton (Goldstone mode) ii) an internal mode which 
appears when the system is perturbed by the external 
magnetic field and iii) the ipk solutions which constitute 
the continuum modes and correspond to magnons. 

In order to find the damping coefficient we must com- 
pute the reflection coefficient H(k). This can be done 
by inserting the numerical results of the even and odd 
phase shifts into the general expression (p8|). In Fig. 3, 
we have plotted TZ(k) for different values of the pertur- 
bation parameter p for the whole range of k. As it can be 
seen the major contribution for the reflection coefficient 
comes from the low energy states, in agreement with the 
well behaved potentials (B5h and (J3q). 



The values of 5 e and S for k — O^are in agreement 
with the ID version of the Levinson'sE3 theorem which 
establishes that 



6«(k = 0) = ttK - ~) 



o 



O " ^0.5 

z a 



S°(k = 0) = 7m°, (37) 

where n e and n° are the number of even and odd parity 
bound states. 



o 



3.0 




' 0.0 1.0 2.0 3.0 4.0 

Momentum 

FIG. 2. The odd phase shift as a function of the momen- 
tum. The continuous line correspond to p = 0.14, the dotted 
line to p — 0.31 and the dshed line to p — 0.50 



' 0.0 1.0 2.0 3.0 

Momentum 

FIG. 3. Reflection coefficient as a function of the momen- 
tum. The continuous line corresponds to p = 0.14, the dotted 
line to p — 0.31 and the dashed line is for p = 0.50. 

Having done that, one can immediately integrate the 
function lZ(k) in expression ( p7j ) which finally allows us 
to describe the damping coefficient as a function of the 
temperature (see Fig. 4). It is important to notice that 
we have not considered the odd bound state of the poten- 
tial ( |34| ) in computing the damping coefficient because in 
evaluating the scattering matrix (Bq), only clastic terms 
are taken into account (see for instance Ref. 9, Ref. 12 
or Ref. 13). 

As it can be seen, the damping coefficient is linear for 
high temperatures. This result can be obtained directly 
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from (^7|). I n fact, for T high enough the damping con- 
stant can be approximated by 



7CO 



1 



2nM s {i J Q 



at, — — — oc 1 

E 



(38) 



which is linear in T, independently of the explicit form 
of 1Z(E). In the low temperature regime we can write 



7(T) 



1 



2vrM, 



dE K(E)(3Ee 



-0E 



(39) 



where E always presents a gap determined by the pres- 
ence of the magnetic field and/or the spontaneous stag- 
gered magnetization. Here we shall not attempt to write 
an approximate expression for (f?9]) because the correct 
behavior of the reflection coefficient was only numerically 
determined. As it is shown in Fig. 4, for low enough tem- 
peratures, the damping coefficient drops exponentially to 
zero due to the existence of the gap. As the temperature 
increases the damping coefficient rises following a power 
law behavior until it becomes linear for high enough tem- 
peratures. This strong temperature dependence of the 
damping parameter, for T below the transition temper- 
ature, will influence directly the correlation function be- 
tween the magnetic solitons. 



B. Soliton mobility for T > T N 

To perform the calculation of the 7r-soliton mobility 
for T > Tn we simply set to zero the b± in the hamil- 
tonian (^J) and therefore, the equation of motion for the 
^-dependent part of the spin degree of freedom (0) be- 
comes a SG equation with solitonic solution ( |l3| ) . In this 
case the potential involved in the Schrodinger-like equa- 
tion ( p7| ) which determines the fluctuations around the 
soliton solution have the form 



o 
U 

C 

■a o.4 
s 



U"{z) = £ 2 (l-2sech 2 £z), 



(40) 



where £ = 6/2. The spectrum of ( fib] ) contains a bound 
state with zero energy 



tpo = A /-sech(£z), 



k 2 



0, 



(41) 



which constitutes the translation mode of the soliton 
(Goldstone mode), and— a continuum of quasiparticles 
modes (magnons) givenli2 by 



k n + it- tanh(£z) 



where 



Itvk S(k n ) 



L 



L 



S(k) — arctan 



2£fc 

2 _ £2 



(42) 



(43) 



As it was already mentioned, the reflection coefficient 
1Z for a general symmetric potential can be expressed in 
terms of the corresponding even and odd phase shifts by 
the relation (S). Re-expressing ( fiHf ) in terms of parity 
eigenstates it is easy to prove that the potential (40) be- 
longs to the class of reflectionless potentials because its 
phase shifts are given by 



S e '°(k) = arctan(£/fc), 



(44) 



that do not distinguish between odd and even parities. 
Therefore no matter how high the temperature rises 
above Tn the damping coefficient is always zero and as 
a direct consequence the ballistic regime for the soliton 
results. 

As we have seen, the solitonic solutions in TMMC 
have different regimes for T below and above Tm- Below 
the transition temperature the 27r-solitons behave like a 
Brownian particle with a finite damping parameter. On 
the other hand, for T above Tjv the 7r-solitons have infi- 
nite mobility corresponding to the ballistic regime. The 
next section is devoted to studying the influence of the 
changes of the solitonic solutions mobility in the dynam- 
ical properties of TMMC. 



IV. DYNAMICAL STRUCTURE FACTOR 



Temperature kT 

FIG. 4. The damping coefficient as a function of the tem- 
perature for different values of p. The continuous line for 
p = 0.14, the dashed line for p = 0.31 and the dotted line for 
p = 0.50 



In this section we will investigate the dependence of 
the dynamical properties of TMMC with respect to the 
temperature and the magnetic field. For T below the 
transition temperature, this will be done through the 
computation of the dynamical structure factor of a di- 
lute gas of 27r-solitons in a dissipative regime. With this 
result, we can analyze the main differences n with the as- 
sumption of ballistic regime using by Holystu in the same 
situation. 
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In a general form, the longitudinal and transverse dy- 
namical structure factors with respect to the external 
field B can be defined as 



(S II(X) = 



f 



{2-Kf 



dtdze l{qz - ut) (S x(y) (0, 0)S X{V] (z, t)) , 

(45) 



where S x ^ corresponds to the spin component in the 
x(y) direction. 

To begin with, let us recall the main results for the lon- 
gitudinal dynamical structure factor reported in Ref. 3. 
Using the model of non-interacting 27r-soliton gas in the 
ballistic regime 5" (q, uj) can be written approximately as 



Sll(^Hn 2 ^ 2 |^(g)| 2 ^, 



where 



0E 2 k [3E 27I uj 2 
exp - 



2ttc 2 



2cV 



(46) 



(47) 



E 27I = 2Bgfi B S 



^l + 4b ± /b 2 

1/-H.-V2, 



-Ab±b~ 2 sinh" i (fe6 ± 1/ 72) 



and 



*£r(?) 



sin(g<i 7I .a) 



2 cosh(g7rd 7r -y/l — a/8) ' 



(48) 



(49) 



a = - In - - 1 + -VT^a 

8 V a a 



d« = -. (50) 
o 



The correlations described by (|46j) are induced by sin- 
gle kinks moving from the origin to the position z in a 
time interval to t. As expected, the Maxellian velocity 
distribution used to describe the 27r-kink gas is directly 
reflected in the Gaussian dependence of the longitudinal 
structure factor with the frequency 

On the other hand, to get a better idea of the changes 
in the dynamical properties when we cross the transition 
temperature, it is convenient to calculate the dynamical 
structure factor for T above the transition temperature. 
As it was demostrated before, above T/v, the 7r-solitons 
moves without dissipation and, therefore, the ballistic 
regime is valid. Using again the model of a dilute gas of 
solitons, the dynamical structure factor can be written 
as 



where 



S' 2 



(2tt) 3 / 2 I/ c 2 q 



§^|i%)| 2 exp(- 



2c 2 q 2 



(51) 



2tt 



F^(q) = — sech( 9 7r/6) and E n = Bg/isS. (52) 



As it can be seen, the dependence with the frequency 
remains almost unchanged no matter what the tempera- 
ture is. At the same time, it should be noticed that once 
the intensity of the central peak and the density of kinks 
in ([l7|) and (|5l]) are proportional, the intensity of the cen- 
tral peak for T < Tjv will be lower. This is a consequence 
of the smaller number of solitons for temperatures below 
the transition temperature. □ 

With the previous results for T above and below T/y 
in mind, we can go further on and study the influence 
of the dissipative regime in the dynamical properties of 
the 27r-kink gas. As it was shown before, below the tran- 
sition temperature the 27r-solitons move in a dissipative 
regime. Therefore, the position of the center of mass for 
those kind of excitations as a function of time can be 
written as 



= 20 + ^7(1- exp - 7 (T)i), 
IK 1 ) 



(53) 



where zq is the initial position, vq is the initial velocity 
and 7(T) is the temperature-dependent damping param- 
eter. Now, to calculate the dynamical structure factor we 
will use the 27r-soliton solutions (|9|)-(|Tl]) with z(t) given 
by ( p3[). Following the same procedure that led to equa- 
tion (|46|) and, after some calculations, iS"(q,w) can be 
written as 



5ll(?,w) 



2n 27r S 2 



\F x Jq)\ 2 T(q,L0) 



(54) 



where 



r = 



E 

n=0 



(-ir 



2f3E 2vl 2 



n 2n 



2m7 



m 2 "/ 2 + lu' 1 



m=0 



and 



C 2 n 
rv, 



(2n) 



m\(2n — m)\ 



(55) 



(56) 



As it can be seen, the dissipative regime for the mag- 
netic solitons in the case in which T < Tjy, changes 
considerably the behavior of the longitudinal dynamical 
structure factor. Although the expression (|54|) is valid 
for all finite values of q, we couldn't perform the entire 
sum to get a closed expression. Therefore, it is helpful 
to study the behavior of (|54])-(|5^) for small momentum 
in order to compare it with the ballistic behavior result 
Assuming that 



g« ^^2/3E 27[ , 



(57) 



the dynamical structure factor S"(q,u>) can be written 

as 



2a 



S\\(q,u,) = —\F^(q)\ 2 A(q,uj) 



(58) 



7 



where 

A(q, u>) — 2it5(oj) exp - 



2 2 



2 2 

g c 

7 



47 2 + 



(59) 



Within the approximation of small momentum, the 
behavior of S^(q,uj) with the frequency, changes from 
the 'Gaussian' central peak to a 'Lorentzian' dependence. 
Therefore, as the temperature is lowered below T/v, the 
central peak behavior is replaced by a smoother-one in 
the frequency domain. This result is a direct consequence 
of the dissipative regime of the 27r-soliton and, as the 
damping constant 7 can be controlled by changing the 
temperature and the magnetic field, a possible indica- 
tion of a non-ballistic regime has been found. Another 
quantity that can be computed in order to get a better 
idea of the influence of the dissipative motion of magnetic 
solitons in TMMC is the T x time of NMR. This problem 
is currently being investigated by one of us. 



V. CONCLUSIONS 



In this paper we have analyzed the possibility of iden- 
tifying two different regimes of motion for the magnetic 
solitons in the TMMC antiferromagnet. We were able 
to show that above the transition temperature T/v, the 
7r-soliton moves without dissipation, even from the field 
theoretical point of view. This result is in complete agree- 
ment with the ballistic regime adopted to understand the 
experimental data reported in Ref. 11. 

On the other hand, for T below the transition temper- 
ature the 27r-solitons in the system move with finite mo- 
bility. Therefore, a dissipative equation of motion has to 
be used in the description of the soliton's center of mass 
motion. This difference in the regime of motion is directly 
reflected in the longitudinal structure factor and, there- 
fore, can be used as an indication of the finite mobility of 
the solitonic solutions below the transition temperature. 
The results presented here could be directly compared to 
the experimental data one may obtain when testing the 
TMMC antiferromagnet in this temperature regime. 

Although the formulation used to compute the damp- 
ing parameter is valid for all values of the external mag- 
netic field, we have restricted ourselves to the study of 
very weak fields. However, our formulation can be used to 
study situations with arbitrarily stronger magnetic fields 
and to other magnetic materials that support solitonic 
solutions without any major qualitative difference. For 
instance, we could treat systems modelled by the 1-D 



Dzyaloshinski-Moriya antiferromagnet which can natu- 
rally be described by a 2DSGO hamiltonian indepen- 
dently of the temperature. 
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